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Abstract. 

The motion of unstable fluid interface due to Richtmyer - Meshkov (RM) instability 
incorporating with density variation has been studied in a spherical target using 
Lagrangian formulation. During the compression in Inertial Confinement Fusion (ICF) 
process, the density of deuterium - tritium (DT) fuel increases 1000 times greater than 
the density of gaseous DT fuel within the core of spherical target. We have extended 
the feature of density variation [PRA,84-Mikaelian & Lindl] in spherical geometry. 
Due to convergent shock impingement, the perturbed interface will be nonspherical 
which leads to the density variation in both radial as well as in polar angle. We 
have shown that the interface of perturbed surface decreases with time to reach a 



Richtmyer - Meshkov instability in a spherical target with density variation 2 
minimum and then kick back to gradual increase. As the perturbed radius decreases, 
the density increases and reaches a maxima corresponding to a minima of perturbed 
radius. This is the practical situation of density characteristics during implosion of 
ICF. The numerical results based on our analytical work show a good qualitative 
agreement with some experimental as well as simulation results. 



1. Introduction 

The direct driven Inertial Confinement Fusion (ICF) process can be divided into four 
stages, i.e., the laser radiation, the implosion-compression, the fusion ignition and the 
deuterium-tritium (DT) fuel burn. Mainly, in the second stage, the hydro dynamical 
interfacial fluid instabilities can act. The Rayleigh- Taylor (RT) instability can occur 
during the two phases of the ablative acceleration and the compression deceleration of 
the ICF implosion. In addition to RT instability, RM instability can occur when shocks 
pass through the irregular fluid interface. The nonspherical implosion in ICF process 
generates the shear flow causing the Kelvin- Helmholtz (KH) instability. The growth of 
these instability can be affected by the density ratio of two fluids, geometry 
of the interface, compressibility effects, heat conduction and mass flow at 
the interface. These instabilities should be mitigated for efficient compression in ICF 
experiment. 
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Mikaelian and LindlJT] explained that the smooth density gradient at the interface 

stabilize the hydrodynamic instability in ICF situation. However, they considered a 

constant and exponential density gradient at the interface. Initially this work started 

by Plesset [2] for stability condition of underwater explosion considering the velocity 

along the radial direction only. He discussed the stability and instability conditions 

depending on the pressure difference at the interface of bubble. Growth of perturbed 

interface due to density [3] — [3] variation has been discussed in different geometry [5] . RT 

growth stabilization effect has been showed very well by hot spot dynamics [5] — [HI] in 

ablation situation of ICF. However, Ramshaw [TT] extended Mikaelian's[T2], [32] work 

using Lagrangian formulation in hydrodynamical model in a spherical geometry. In 

this context, several authors [Hj — [18] pointed out that during implosion(acceleration 

phase) the fluid is compressed hence the fluid density increases rapidly inside the 

sphere. Similarly, during the blow off(deceleration phase) the fluid is relaxed to expand 

and fluid density falls exponentially from the perturbation surface. During robustness 

compression [HI] — [ZE] of direct driven ICF target, the density can be compressed by 

100 to 200 times of liquid DT density at cryogenic temperature [22] — [23]. Perturbation 

amplitudes and its growth rates of hydrodynamical instabilities of laser fusion targets 

has been studied in acceleration and deceleration phase including the heat conduction 

effect [23] — [25]using Lagrangian formulation. Numerical, simulation and experimental 

results [25] — [2H] °f RMI are performed in cylindrical geometry considering implosion 

and explosion situation. However, the growth of RM instability can be suppressed due 
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to different shape of the interface [30J — [3 lj . Recently, X-ray and proton radiography 

technique has been used to study the fast ignition implosion in cylindrical geometry 

|32] — [31]. The simulation result well agree with the experimental result. So, density 

variation and sphericity are two key factor to understand clearly for successful of ICF 

process. 

In this paper, we investigate the nonlinear evolution of perturbed interface in 
spherical geometry considering mean interface [R{t)} and angular amplitude both are 
dynamical variables. As a result, the growth and growth rate of nonspherical interface 
due to implosion has been changed for both variable with time. Hence, the nonspherical 
irregular interface changes explicitly with radius, polar angle and time. Analytically and 
numerically we have shown the density variation for compression-blow off situation of 
ICF. Theoretically high compression may be achieved depending on the density variation 
at the interface of two fluids. Polar plot also shows the time evolution of the interface in 
spherical geometry. Growth reduction has been shown for high value of mode number 
for a particular initial conditions. 

The paper is organized as follows: section 2, deals with basic assumptions and 
equation of continuity. Expressions of kinetic energy of two fluids followed by Lagrangian 
equation of motion have been discussed in section 3. The nonlinear equations which 
describe the growth and the growth rate of perturbed interface are presented in section 
4. Section 5 contains results and discussions. Finally, we conclude the results with a 
hope of future study on this work. 
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2. Basic equation 

We have considered two concentric spherical shells having two incompressible, 
immiscible and inviscid fluids of density p\ and p%. The unperturbed interface at any 
time is given by r = R(t). The lighter fluid is bounded in the inner shell at R\ < r < R(t) 
and the heavier fluid is ceiling at radius R2, where, R\ <C R(t) <C R2, as shown in Fig.l. 
Here we assume that the interface is perturbed to r = f which is given by 



r(t,6) = R(t) + J^lla^PiicosO) 



r(t,9) = \R{t) + 8R ai ] + J^la^PiicosO) (1) 



where, ai <^ R and R,ai, Pi(cos6) are mean radius of perturbed interface, angular 
displacement and Legendre polynomial of mode number / respectively and 5R ai = 
R[l — (^£) 2 ]. Here we ignore the azimuthal dependence of perturbed interface due 
to spherical symmetry. However, simulation [35] study shows that there is a weak 
dependence of growth rate on azimuthal angle. 

Now we assume that the fluid motion is irrotational in the interface region,i.e., its 
obey Laplace's equation V 2 $ = 0. So the velocity potentials for lighter and heavier 
fluids are 



*, = + I J?i±l (Rd, + 2ft.,) P,(cvs0) (2) 



respectively. 

Due to shock impingement, the nonspherical perturbed interface leads to the density 
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variation as follows 



pi(r,6,t) = pi + pu(r,t)Pi(cos6); 



Ri <r <r 



(4) 



p 2 (r,9,t) = p 2 + p 2 i(r,t)P l (cos9); 



r <r < R 2 



(5) 



for lighter and heavier fluid, respectively. This type of density variation feature was not 
considered in previous studies. 



2.1. Equation of continuity 



The equation of continuity in spherical polar coordinate can be written as 
dt r dr r dt 



(6) 



where, u ir = ^ and u ie — are radial and transverse component of fluid velocity, 



de 



respectively. Substituting the value of U\ r , w^and p 1 , we get the following equation for 



lighter fluid 



dp 



ii 



dt 



R 2 R 21+1 



(i-i)' 



(Rd l + 2Ra l ) (-^r) Pi(cos 



dp 



ii 



dr 



+ W«+i (ffil+2AlI )f^)§ (CM9) 



I 



o 



(7) 



Neglecting the product term (a; x p n ), we get the following partial differential equation 
dpn R 2 R dp n 



dt 



+ 



r 2 dr 







(8) 



The solution of the above p.d.e is given by, 



Pn(r,f) = piies 



[il 3 (t)-r 



(9) 



Similarly, we get the density variation for heavier fluid as follows 



P2i(r,t) = p 2 ie 3 



f [r 3 w- 



(10) 
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where, p is a positive constant for implosion situation. 

The mass conservation at the interface gives us the following relation 

(hi _ P21 

Pi P2 

3. Kinetic energy 

To study the equation of motion of perturbed interface we have calculated the kinetic 
energy of the two fluids. 
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(11) 



3.1. Kinetic energy of lighter fluid 



The kinetic energy of lighter fluid can be calculated in the following way 



T 1 = ±JdnJ^ Pl (r,6,t) 



|V$i| 



r dr 



Pi + Pu(r,t)Pi(cosi 



1 + ' 1 



dr 



- fdn r 

2 J Jr! 



HLti ] [Rd l + 2Ra l 



— Tio + TirQ 



Pi + Pu(r,t)Pi(cos9) 

\ 2 / r 2(l-l)\ 



R^R 2 R 2 R 21 + 1 



R 21 



+ 2 



Pi(cos6) 



1 (21 + 1 



P 



r 2 \ 06 



Rdi + 2Rai 



Rd t + 2Ra x 



r 



r 2 dr 
(i-iy 



R l 



Pi(cos6) 



2 f r 2(l-l)\ ( d p\ 2 - 



R 21 / V d9 



r 2 dr 
(12) 



Where, T 10 and T lr0 are the kinetic energy coming from homogeneous density p\ and 
inhomogeneous density p u (r,t) part of the lighter fluid, respectively. Now using Eq. 9, 



we get Ti rd as follows, 



Tire = ^f J Pi(cose)<mf 



R 4 R 2 



+ 2R 2 R\ 



'21 + 1 



+ 



Ri 

(21 + 1 



Rdt + 2Ra t 



~RJ~ 



Pi(cos6) 



V 2 



Rd l + 2Ra\ (^\ (Pi(cos6) 



Richtmyer - Meshkov instability in a spherical target with density variation 

ft 3 (t)-r 3 



air (13) 



Similarly, we have calculated the kinetic energy for heavier fluid in the preceding section. 



3.2. Kinetic energy of heavier fluid 



T 2 = ljdnJ% 2 (r,9,t) 



|V$ 2 | 2 



r 2 dr 



P2 +P2i(r, t)Pi(cos6) 



dr 



r 2 \ 86 



r dr 



1 r rR2 
- / dtt 



P2 + P2i(r, t)Pi(cos9) 



R 4 R 2 R 2 R 2/ + l/„. ^ \(R^ l+i y 



+ 2- 



Rdi + 2Ra 



/ \r 



(1+2) 



^4^) ( Rd i + 2Ra i] ( ^LJr ) {Piicost 



.2(1+2) 



21 + 



(J + l) 



r dr 



— T20 + T 2r 6 



Where, T 20 and T 2r9 are the kinetic energy coming from homogeneous density p 2 and 
inhomogeneous density p 2 i ( r > t) part of the lighter fluid, respectively, and / dVt = 



2n J sin 9d9=so\id angle. Now using Eq. 10, we get T 2r e as follows, 

( li 1 ' ■ \ 



'- 2r0 



P21 



/rR-2 



R^R 2 „R 2 R 21 + 1 , _ , , ... , , ., 

+ 2 — — ^ — ( Ra t + 2Ra t ) [ ) P^cosO) 



+ 



r^ r 
(21 + 1 



V 2 
21+ r 



Rdi + 2Ra 



Vi? 2( ' +1)x ' 



r 2(l+2) 



Pi(cos9) 



2i 



[&(t)- 



Jr 2 rfr (15) 



3. 3. Lagrangian equation of motion 



Now we have formulated the problem using Lagrangian equation of motion. The present 



problem is related to the RM instability i.e., shock impingement at the interface of two 
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fluid and hence neglected the potential energy. So, we have Lagrangian equation of 

motion as follows 

where, T = T\ r o + T 2r g and q — > R and a\ both are dynamical vriable. 



4. Nonlinear equation 

Using eqs.13 & 15, We get from eq. 16, a pair of coupled nonlinear second order following 
differential equation. The detailed calculations are given in appendix. 

f 1R (R, ai )R + f 2R (R, a^at + f 3R (R, a^R* + f 4R (R, a^atf + f 5R (R, a^Rdt = (17) 

haXR^a^R + haXR^a^ + f 3ai (R, ai )& + f 4ai (R, a^atf + f 5ai (R, a^Rcn = (18) 

From this two second order nonlinear differential equation we get following four first 
order differential equation, which describes the growth and growth rate of perturbed 
interface. 

alX 1 



= X 3 (19) 
= X 4 (20) 

(21) 
(22) 



dr 

dX 2 
dr 

dX 3 _ n\d 2 — n 2 d\ 
dr m\n 2 — m 2 n\ 

dX A m 2 d\ — m\d 2 



dr m\n 2 — m 2 n\ 

where, X 1 = ^,X 2 = j£ m 1 = £ 1 (X 1 ,X 2 ),m 2 = &(X 1 ,X 2 ),n 1 = ^(X U X 2 ) and n 2 = 
^ 1 (X 1 ,X 2 ) are nondimensional perturbed radial displacement, angular displacement, 
radial velocity and angular velocity respectively. 



Richtmyer - Meshkov instability in a spherical target with density variation 10 
5. Results & Discussions 

The feasible analytical solution of this set of equations is quite impossible. We can 
numerically solve this set of equations using Runge-Kutta-Fehlberg numerical technique. 
The numerical results shows that due to the shock impingement, the interface decreases 
and consequently the DT fuel fuel is compressed during acceleration phase with time. 
The perturbed interface attains a minimum position and then kicks back to gradually 
increasing in time (Fig.2a) during the deceleration phase. The density increases as 
the perturbed interface decreases and attains a maximum density corresponding to 
minimum of perturbed interface then the density gradually decreases (blow-off situation) 
as the perturbed radius increases again(Fig.2e). This is the practical situation of 
compression-blow off process of ICF. During acceleration phase the fuel compression 
is maximum. The dotted line which also represents the growth of interface with time 
when a = nR\ = 0. However, dash-dot (blue) and solid (red) line for a = 10 and 15 
respectively, represents the stabilization due to density variation. Our this analytical 
results show a qualitative good agreement with simulation [7] as well as experimental 
results [32], [S3], [34J. Though the experiment has been done in cylindrical geometry, the 
radius of a cylinder can be equivalently represented by the radius of a sphere. We have 
also shown the polar plot of perturbed surface at different time (Fig. 3, 4, 5). We have 
considered the density ratio r = 10, 1 = 50(Legendre polynomial mode number) for 
numerical calculation. This type of density ratio [13] is appropriate for ICF process. We 
have shown in Fig. 6 that the growth rate is maximum at I = 50. 
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6. Conclusion &z future scope 
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Maximum compression can be achieved depending on the density variation at the 
interface. The acceleration and deceleration phase can be better understood from the 
numerical results. However, in this work we do not consider the temperature effect 
which is more realistic. In future we can take attempt to include the temperature effect 



also. 
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Appendix 



The details of calculation for kinetic energy is given below. 

~R 4 R 2 _ „o ^ 121 + 1 ' 



T ire = ^f J Pi(cos6)dnJ^ 



+ 2R Z R\ 



Ri 
21 + 1 



Rdi + 2Rai 



Pi(cos9) 



2 ^Rcn + 2Ra^ (^j (Pi(cos9)j 



R A (t)-r [i 



dr 



p +(<m" 



Ri 



R 4 R 2 



-Pi + 2R 2 R 



+ 



^ 2l + 1 ^^Rd l + 2Ra^j (^j 



V 2 



121 + 1 
^) 3 + 



Rdi + 2Ra x 



r 



R l 



(Pif 



, 1 21 + 1 
2 I 2 



— Tip! + T\p 2 + Tip 3 + Ti Pf) 



Where, 



R 4 R 2 



ipi 



Pn / Pidn 



r | m 

o 3 



R :i (t)-r 3 



dr 



Ri r z 



(25) 
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. \ /rC- 1 ) 



■ip 2 



2 J Jr 1 



n ■ 21+1 

2R 2 R\ 



Rdi + 2Rai 



R l 



(PiY 



R 3 (t)-r 3 



dr 



Ri 



Ri 



'21 + 1 



2 / \ 21 



Rd l + 2Ra l j (P t f 



R 3 (t)-r 3 



dr 



R 3 (t)-r 



dr 



12 
(26) 

(27) 

(28) 



After a straight forward calculation of Eq. (25) and putting the value of f and expanding 



it upto second order of (^), we get the following equation 



(29) 



Similarly, 



121 + 1 



T 1P2 = AnR 3 Pl ) ) (ifci, + 2Ra t ) 



(Z + 3 + ^i? 3 ) (a,/i2) 



+ 



T 1Ps =4rf ft (^)(?i±i)( ffi ' + 2 ^ 



i(2/ + l)(/ + 3) ,/ 2 (2/ + l) 
A^ 3 \ 



(30) 



fiR\X l 



>/2(2/ + l)' 



2(/ + 2)7 (2/ + l) 2 
1 1 / a A 2 



2 ' 2 (2/ + 1) 2 J \ yJ R / 



(31) 



T lPe =4tt j R 3 pi(— )(2/ + l 



Pi 



2lR 



/(/ + 1) 



+ 2 + ^)r, 



(2/ + l) 2 (/ + 2) 



^(>. + {^('-^-^}'('-)(r 



P21 



2rt9 



/[■Ri 



Rdi + 2Ra t 



2(2/ + 1) 



(32) 



(21 + 1 



Rdi + 2Ra x 



Pi 



1 (2l+l\(„. . „a VfR 2 ^ l+r >\(dPi\ 2 ' 



1 I r 2(l+2) 

^2Pi + T 2 p 2 + T 2 p 3 + T 2 p e 



e^~ r % 2 dr (33) 
(34) 
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Where, 



t 2P , = ( f i (i + ^y^L Ql - + 



2 r 



fiR\( ai /R) 2 ^ ( ai /R) 



v / 20 + i) 
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Figure 1. Geometry of the problem 
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Figure 2. 
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Polar plot of perturbed radius at x=4 
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Figure 3. 
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Polar plot of perturbed radius at x=6 
Red and blue line for a=10 and respectively 
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Figure 4. 




Figure 5. 
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Mode no.(l) vs. Growth rate plot 
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Figure 6. 



